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ABSTRACT 

We use the non-linear spherical model in cold dark matter (CDM) cosmologies with 
dark energy to investigate the effects of dark energy on the growth of structure and 
the formation of virialised structures. We consider dark energy models with a constant 
equation of state parameter w. For —1 < w < —1/3, clusters form earlier and are 
more concentrated in quintessence than in ACDM models, but they form later and 
are less concentrated than in the corresponding open model with the same matter 
density and no dark energy. We point out some confusion in the literature around the 
expression of the collapse factor (ratio of the radius of the sphere at virialisation to 
that at turn-around) derived from the virial theorem. We use the Sheth & Tormen 
extension of the Press-Schechter framework to calculate the evolution of the cluster 
abundance in different models and show the sensitivity of the cluster abundance to 
both the amplitude of the mass fluctuations, as, and the as — w normalisation, selected 
to match either the cosmic microwave background observations or the abundance of 
X-ray clusters. 

Key words: cosmology: theory - galaxies: clusters: general - large-scale structure of 
universe - galaxies: formation 



The cosmological picture that has emerged from a host of 
recent observations is that of a spatially flat universe domi- 
nated by an exotic component with positive energy density 
and negative pressure, the dark energy, which is responsible 
for the observed accelerated cosmic expansion. 

The nature of dark energy remains unknown. It 
could be Einstein's cosmological constant or the man- 
ifestation of a scalar field slowly rolling down its po- 



1988, 



1998, 



tential, as in quintess e nce models iPeebles fc Ra~ 
iRatra fc Peebles! 11983. ICaldwell. Dave fc Steinhardt 
see iPeebles fc R atra 2003 for a review). The equation of 
state, w — px/px, where px is the pressure and px the 
density of the dark energy (hereafter referred to as the 
X-component), describes how the X-component evolves as 
the universe expands. Determining the equation of state 
will help reveal the nature of dark energy and discrimi- 
nate among models, as the cosmological constant is char- 
acterised by a constant w = — 1 whereas in quintessence 
models w can be different from —1 and vary with cosmic 
time, w is directly related to the shape of the potential of 
the quintessence field. It is a fundamental parameter that 
remains poorly constrained observationally, partly because 
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of degeneracies with other parameters, mainly fio and the 
r.m.s . of the amplitude of fluctuatio ns within 8 h~ Mpc, as 
(e.g.. lLokas. Bode fc Hoffmarfeooj) . 

Dark energy not only affects the expansion rate of the 
background and the distance-redshift relation, but also the 
growth of structure. The collapse of overdense regions due 
to gravitational instability is slowed down by the Hubble 
drag due to the expansion. This effective friction depends 
on the expansion of the background. The formation rate 
of haloes, their evolution and their final characteristics are 
modified. Dark energy is therefore expected to have an im- 
pact on observables such as cluster number counts and lens- 
ing statistics due to intervening concentrations of mass on 
the line-of-sight of background sources. 

The abundance of rich clusters of galaxies can be 
used to constrain the cosmological model and the prop- 
erties of dark energy. This has been discussed first 
for pure cosmological constant models fe.g.. [Lahav et alJ 
I l99ll lLilidll99l lLacev fc Coldll99l IViana fc Liddldll996L 
lEke. Cole fc Frenkl 11996) and then for models with dark 
energ y where w could be differe nt from —1. In a seminal 
paper IWang fc Steinhardt] lll998l) have discussed the evolu- 
tion of the cluster abundance in general quintessence models 
with n e gative pre s sure. Following up on that paper, lLokasI 
1200 ll) . Basilakosl £2003) have examined cluster formation 



in models with constant w. iMota fc van de BruclJ (1200 
have selected par t icular potentials for the quintessence field. 
iBattve fc Welled (2003) have investigated the constraints 
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on dark energy from future Sunyaev-Zeldovich surveys. 
lLokas. Bode fc Hoffman! i2004h have used JV-body simula- 
tions in flat models with constant w to measure cluster mass 
functions. 

In this paper, we first discuss the linear evolution of per- 
turbations in dark energy models with constant w. Then, we 
use the non-l inear spherical mod el ("top-hat" model first 
developed bv iGunn fc Gott lll972t). Rather than using the 
approximation given by Wan g fc Steinh ardt ( 1998) which is 
strictly valid for w = —1, we calculate explicitely the po- 
tential energy associated with the dark energy component 
inside the collapsing sphere and the ratio of the radius of the 
sphere at virialisation to the turn-around radius. We follow 
the evolution of overdense regions as they first expand with 
the background before collapsing and settling into virialised 
structures. We then calculate t he ab undance of massive clus- 
ters in the|Preg s & Schcchtcr ( 1974) framework. We use the 
ISheth fc Tormerall999r) extension which has been shown to 
be in remarkable agreement with results from A^-body sim- 
ulations of hierarchical structure formation. 

We use five cosmological models, all with Ho = 70 
kms _1 Mpc~ 1 : three dark energy models with the same en- 
ergy content: fio=0.3, fix,o=0.7, but that differ in their 
value of w: ACDM with w — — 1, a w — —0.8 model, a 
w — —0.6 model, and for comparison an Einstein-de Sitter 
model with fio = 1 and an open model with fio = 0.3. 
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Figure 1. Evolution of the scale factor in the different models. 
In quintessence models, the evolution of the scale factor is inter- 
mediate between that in a ACDM model (continuous line) and 
that in an open model (dashed-dotted-dotted-dotted line). 



2 THE BACKGROUND COSMOLOGY 



We consider a cosmological model with two components: (1) 
a non-relativistic component which comprises all forms of 
matter (luminous and dark) that clusters under the action 
of gravity, and (2) a uniform X-component with negative 
pressure which does not cluster at the scales of interest. In 
general, quintessence models have a component that is spa- 
tially inhomogeneous but does not cluster on scales less than 
100 Mpc. The Friedmann equation that describes the evo- 
lution of the scale factor a with cosmic time is modified to 
include the effect of the X-component: 



H(a) 
Ho 



O P ^ J-O n- 2 MCI PX ^ 

ilo h sifcffl + «x,o- 



Px,o 



(1) 



where H(a) — a/a is the Hubble parameter, fi are the den- 
sity parameters (ratios of the energy density to the crit- 
ical density p cr it = §|^j) for the matter, the curvature 
(flfc = 1 — fio — fix.o) an d the X-component respectively. 
The parameters are indexed by when they refer to the 
present time. 

The energy density of the X-component varies as 

px,o 

with n = 3(1 + w). 

Different n and w correspond to different types of en- 
ergy densities, with n = (w = — 1) corresponding to the 
cosmological constant. Note that n doesn't need to be either 
integer or constant. It may vary with cosmic time, as in ki- 
nessence models. In the following we consider only models 
with constant n and w. 

Note also that for n = 2 (w = —1/3), the dark en- 
ergy term can be combined with the curvature term in the 



Friedmann equation; flx,o cancels out and the expansion is 
the same as for a flo universe. Thus a flat universe with 
w = —1/3 behaves like an open universe without dark en- 
ergy and with the same Qq. Cosmic acceleration (i.e., the 
existence of an inflexion point in the a(t) curve) requires 
n < 2. 

The Friedmann equation can be written more simply as 



Ho 

a = W) 

where 



(3) 



f(a) = a 1 [Q. a 3 + (1 - Q - Qx,a)a 2 + Qx,oa "] 



-1/2 



(4) 



Integrating the Friedmann equation, one obtains the 
variations of the scale factor with cosmic time in different 
cosmologies (see Fig. . Universes with a dark energy com- 
ponent are older than with fix = and give perturbations 
more time to grow. For the currently favoured parameters 
the evolution of the scale factor in quintessence models is 
intermediate between that of a cosmological constant model 
(w = — 1) and of an open universe. 

The density parameters vary with redshift z = l/o — 1 



Qo(l + zf n x , {l + z) n 
"(z) = — t^-, — ; yi x {z) = 7$T7-\ (5) 



EHz) 



E\z) 



where E(z) = Z&L. In the case of the cosmological constant 
commonly denoted A, 



fix = A 



3H 2 



(6) 
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Figure 2. Evolution of the growth factor normalised to the 
present time and divided by the scale factor a as a function of 
redshift for different cosmological models. In an Einstein-de Sit- 
ter universe -D+/a = 1 at all redshifts. D+/a is highest in the 
open model. In quintessence models it takes intermediate values 
between the open model and the ACDM model. 



3 GROWTH OF LINEAR PERTURBATIONS 

In linear theory, the evolution of the density contrast of a 
perturbation 5 = (p pe rt — p)/p is governed by 



5 + 2-5- 4nGpS = 
a 



(7) 



which can be solved numerically. Anal ytical s o lution s 
have been given for Ap = (e.g. iPeeblesI ll98Cft. 
Op + Ap = 1 ij Bikiliauer^Bucher^^^^asail ^9^: 
[Chernin. Nagirner fc StarikovaT 120031: lLokas fc Hoffman! 
2001). A way to estimate the growth f actor numericall y 
for any (Op, Ap) pair has been given by lHamilto'nl i200ll) . 
iBasilakoa i2003h provides a generic approximation formula 
for D+ (z) as a function of Sip , Qx and w that agrees 
with the approximations g i ven b y lLahav et alJ dl99ll) 
and ICarroll. Press fc Turner! lll992t) for the ACDM case. 
iLinder fc^eid{m?T2"o"o"3Tl discuss the influence of the equa- 
tion of state (including a time-varying w) on the growth of 
structure. 

It is common practice to normalise the growth factor to 
its value at the present time. Fig. [5] shows the variations of 
the growth factor normalised at z — and divided by a as 
a function of redshift for the different cosmological models. 
In an Einstein-de Sitter universe D+/a is always equal to 
1, and structures can form at any redshift. The curves show 
that the growth factor is highest in an open universe, and it 
takes intermediate values for decreasing w up to w — — 1, the 
cosmological constant case. A ACDM universe has a higher 
growth rate than an Einstein-de Sitter universe throughout 
its history, but the growth rate drops at redshifts below ~ 
1/fio- For —1 < w < —1/3 the decrease in the growth rate 
occurs at higher and higher redshifts and over a large range 
of redshifts as w increases. In other words, structure growth 



will be stronger at early times in an open universe compared 
to a ACDM universe, with intermediate rates in quintessence 
models. But structure growth is reduced earlier in open and 
quintessence models, although it is still significantly higher 
than in ACDM models. 



4 NON-LINEAR EVOLUTION OF A 
SPHERICAL OVERDENSITY 

The equation of motion of a spherical shell in the presence 
of dark energy is: 

r 4:nG 
r = 3~ 



(/Ocluster + PX.eff) 



(8) 



where p c iustcr is the time-varying density inside the forming 
cluster, and 



px,eff = Px + 3p x = {n — 2)p x 



(9) 



is the effective energy density of the X-component. px.eff is 
constant for the cosmological constant; in a universe with 
a dark energy component with n — 2 (w — —1/3), pertur- 
bations evolve in the same way as in a universe with the 
same f2p and no dark energy. Note that for n^0 and n/2 
(w ^ — 1 and w 7^ —1/3), the density of the X-component 
inside the overdensity patch is dependent on the evolution 
of the background. This makes it impossible to integrate di- 
rectly the equation of motion and cast it into a first-order 
energy equation. The second-order equation of motion has 
to be integrated taking into account the time variation of 
the dark energy. 

We shall follow the evolution of a spherical overdense re- 
gion with some initial overdensity. At early times, it expands 
along with the Hubble flow and density perturbations grow 
proportionally to the scale factor. If the initial overdensity 
exceeds a critical value, the overdense region will break away 
from the general expansion and go through three phases: 

(i) expansion up to a maximum radius; 

(ii) collapse; 

(iii) virialisation. 

We call Zta the redshift at maximum expansion (turn- 
around) and Zcoii the redshift at which the sphere suddenly 
virialises, r ta and r v i r the corresponding radii of the sphere. 
Ideally, the sphere should collapse down to an infinitely small 
radius, but it is assumed that this is prevented by the growth 
of small inhomogeneities. Let us discuss in turn the different 
phases of evolution of a collapsing sphere. 



4.1 Expansion up to rta 

Let us calculate the overdensity of the forming cluster at 
turn-around 



c = 



pclustcr 



Pbg 

Let us define 



a 

Ota 



and 



y 



r 



(10) 



(11) 



where ata is the scale factor of the background when the 
perturbation reaches turn-around. For a flat universe, the 
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Figure 3. Overdensity at turn-around versus collapse redshift in 
different cosmological models. Clusters are denser at turn-around 
in quintessence models than in the ACDM model. Quintessence 
models are intermediate between the ACDM model and the open 
model. Q ~ 5.6 in an Einstein-de Sitter universe. 



evolution of the background and that of the perturbation 
are governed by the two following equations: 



x = H t ^li 2 [n(x)x]- 1/2 



7 + <- 2,! §r* 



(12) 
(13) 



where we have used the fact that the mass of the forming 
cluster is conserved: 



3 _ 3 

Pclustcr^ — Pclustcr,ta^ta • 



(14) 



Note that the second term in eq. 1131 vanishes in two cases: 
for perturbations in a n = 2 universe and for the ones in 
a universe with fix ~ 0. In the n = (the cosmological 
constant) case also, the evolution of the perturbation is in- 
dependent of that of the background. 

£ can be determined by integrating the above differen- 
tial equations using the boundary conditions (dy/dt) x =i = 
and (y) x —o = 0. The variations of C with collap se redshift 
are shown in Fig. [3] W ang fc SteinhardU (|l998) provide a 
fitting formula for £ as a function of f2 ta and w for a spa- 
tially flat Universe. Perturbations collapsing at a certain red- 
shift are denser at turn-around relative to the background 
in a universe with quintessence than in a ACDM model. At 
high redshift, ( tends toward the fiducial value of 5.6 for 
an Einstein-de Sitter universe (as fi tends toward 1 and fix 
toward 0). 

Fig. [I] shows the evolution of a perturbation collaps- 
ing now (zcoii = 0) in different cosmological models. Per- 
turbations reach turn-around and collapse earlier in the 
quintessence models than in the ACDM model, and even 
earlier in the OCDM model. 
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Figure 4. Evolution of a perturbation collapsing at the present 
time (z co n = 0) in different cosmologies. The radius r has been 
normalised to the turn-around radius rt a . The x-axis shows 
the time normalised to the turn-around time for the ACDM 
model. Perturbations reach turn-around and collapse earlier in 
the quintessence model than in the ACDM model and even ear- 
lier in the OCDM model. 

4.2 Collapse and virialisation 

The virial theorem can be expressed in the general form 
T = pU /2 for a system with a potential energy of the form 
U oc R p and kinetic energy T jLandau fc Lifshitzlll96fl(l . For 
the cosmological constant and the quintessence, the poten- 
tial energy is proportional to R 2 , and we can therefore use 
the virial theorem. Combining the virial theorem and the 
conservation of energy one obtains a relation between the 
potential energies of the collapsing sphere at turn-around 
and at the time it virialises: 

Ugm + Uxm = i(7G,vir + 2[/ X ,vir (15) 

where Uq is the potential energy of the matter (Uq = 
— | GI £ ) and Ux is the potential energy of the quintessence 
inside the collapsing sphere. 

The potential energy associated with the quintessence 
Ux can be obtained from the Poisson equation with the 
pressure term: 

A$x =4nG{ Px + 3px) (16) 
which yields: $ x = (n - 2)p x 2j ^-R 2 and 

Ux = {n-2) PX ^R 2 . (17) 

Replacing the energies by their expressions in Eq. (1151 . we 
can obtain a relation between the collapse factor, defined as 



V = 



'ta 



(18) 



and the magn itude of the da rk energy at turn-around. Fol- 
lowing [ilfeT^^hj^ir^ (l200ll) . let us define a quantity 6 to 
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Figure 5. Ratio of the density of the virialised sphere to the 
density at turn-around p v ir/pta = V~ 3 as a function of 6, the 
strength of the dark energy at turn-around, for different models. 
The dotted line corresponds to an Einstein-de Sitter model (den- 
sity ratio of 8). Values are shown for increasing aco11 from 1.60 to 
1.70 by steps of 0.02. 



describe the importance of the X-component relative to the 
gravitational attraction of the matter at maximal expansion: 



/ PX,oS \ 

\ 2p J 



2P 

where px,efi is defined in eq. 
Noting that 



PX.cff 



2 ^3 ^OcollV 

V a ta / 



we obtain the following relation between 77 and 9: 



AO 



( fflcoll 

V ata 



2(1 + 0)77 + 1 = 0. 



(19) 



(20) 



(21) 



We would like to point out that there is some confu- 
sion in the literature around the exac t expression of the cu- 
bic equation H21II. Som e authors (e.g.. llliev fc Shapiroll200ll 
iBattve fe Welleil l200St) have derived the same expression 
as ours, also using the Poisson equation with the pressure 
term f or the dark energy. O t her authors, however, (e.g., 
Lokas] I2001L iBasilakosI 120031 . IWeinberg fc Kamionkowskil 
2003) seem to have used the a pproximation of r v i r /r t a given 
bv lWang fc Steinhardd Jl998t) , which is exact for the cosmo- 



logical constant case, but differs significantly from the an- 
alytical expression when w ^ —1. Rewriting equation (12111 
using the same notations as lWang fc Steinh ardt (1998J), we 
have: 



2 VvV s -2(1 + ^)77 + 1 = 
with 



(22) 



Vt=- 



(n-2) Ox.ta 
(n-2) flx.vir (l + 2 co „) 3 



Note that those eq uations differ from those of 
IWang fc Steinhardtl lll99Sft where the term -(n-2) in the 
expressions above of r/t and rj v is replaced by 2. The formulae 
are identical for n = (the cosmological constant case) but 
they differ when This results in different values of the 

density of the virialised clusters , as we shall see. In partic - 
ular, using the approximation of IWang fc Steinhardd (Il998tl 
in the w = —1/3 case yields to erroneous results. Because of 
the contribution of the pressure term in the Poisson equa- 
tion, a factor (n — 2) appears in the cubic equation, which 
means that for n = 2 (w = —1/3) one should recover the 
same expression for density contrast at virialisation as for 
an open universe with the same fio- This is not the case if 
the (n — 2) term is omitted. 



4-2.1 The collapse factor versus 9 and 

Let us solve eq. 121H to express the collapse factor 77. For 
9 — (which corresponds to px = or to n = 2), 77 = 
1/2. This is a well-known result for models without dark 
energy: the final radius of a virialised sphere is half the turn- 
around radius. For 9 7^ one obtains a cubic equation in 
77 that can be solved analytically to express the collapse 
factor as a function of 9 and for different values of n. 
Here the discriminant of the cubic equation is negative and 
there are three unequal real solutions. The first solution is 
always greater than 1 for the range of 9 of interest, while 
the second solution is negative. The third solution tends to 
1/2 for 9 — > and it can be expressed as: 



77 = 2 
with 



60 



'/acoii\ n/2 / 

-UrJ cos l 



■ + 4n 



(69) 



3/2 



(1 +6») 3 / 2 



V a ta / 



1/2 



(23) 



(24) 



C 



(l + Zta) 3 



For Txi = —1 (the cosmological constant), 77 is independent 
of %^ and i t can b e calculated readily, as was first done by 
lLahav et ail dl99ll) . who also gave a simple approximated 
solution valid for a first-order expansion of rvir/na around 
1/2 (the value for A = 0). 

In quintessence models, 77 depends on ^ aSL , which will 
be calculated below. Fig. shows the ratio of the density of 
the virialised sphere to the density at maximum expansion 
versus 9, the strength of the dark energy at maximum ex- 
pansion, for different cosmological models. In order to reach 
an equilibrium in the presence of dark energy, the sphere has 
to achieve a higher density as 9 increases. The final density 
of the virialised object is even higher in quintessence models 
than in ACDM models, and it increases with increasing w. 
The curves corresponding to quintessence models are inter- 
mediate between those of the ACDM model and the open 
model. 

The minimum radius of the virialised sphere is reached 
for 9 = 0.5. In a universe with a cosmological constant, the 
minimum radius is 77 m i n = 0.366. In quintessence models, 
the radius of the virialised sphere tends toward a minimum 
7/min = 0.333 as the ratio ^ — increases. The maximum ratio 
of the density of the virialised sphere to the density at turn- 
around is thus 8 for models without dark energy, 20 for A 
models, and 27 in quintessence models (when w — ► —1/3). 
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Figure 6. Overdensity A v i r versus collapse redshift for differ- 
ent cosmological models. Quintessence models are intermediate 
between the ACDM model and the open model. 



Expanding r\ around 1/2, its value for 9 — 0, one obtains 
a first-order approximation: 



1 



T) = 



•*te)" 



'(ate! 



(25) 



4.2.2 The ra tio 

The ratio can be derived from the integration of the 
Friedmann equation (J^J to calculate the scale factor as a 
function of time, using the fact that the collapse time is 
twice the turn-around time: 

tcoll _ g 

Oooil and a ta are related by: 



(26) 



f(a)da 



f(a)da. 



(27) 



Although is constant for an Einstein-de Sitter uni- 
verse (equal to (hsll) 2 / 3 — 2 2/3 = 1.587), it varies in other 
cosmologies, decreasing with increasing z co ii, and converges 
toward the Einstein-de Sitter value at high redshifts. 



4.3 Density contrast at virialisation 

In Sect. l4.ll we have calculated the overdensity of the forming 
cluster at turn-around: £ = (pciuster/pb g ) ta - In Sect. 14.21 we 
have used the virial theorem to express p v ir/p c i uster,ta cLS a 
function of (px/pciuster)ta- We are now able to calculate the 
density contrast of the virialised cluster as a function of the 
collapse redshift: 



Avir(Zcoll) 



Pvir 
Pbg.vi] 



V 1 + 2coll / 



(28) 



Figure 7. Decrease of erg with increasing w, using the X-ray 
cluster normalisation or the CMB normalisation. 



fsee | La cey fc Cold |l993l for a A = universe, 
lEke. Cole fc Frenkll996] for a flat universe with a cosmolog- 
ical cons tant. lLo ka*H%^Ioffmanl l200ll for a non-flat universe 
with A, IWang fc Steinhardd Il998l for a flat universe with 
quintessence). The result is shown in Fig. El in quintessence 
models, the overdensity at virialisation takes intermedi- 
ate values between the ACDM curve and the open model 
curve (which is the limiti ng case when w — > —1/3). The 
IWang fc Steinhardtl lll998l) approximation yields an overes- 
timate of A v j r for w > — 1 by up to 40% when w — > —1/3. 
IWeinberg fc Kamionkowskil {2003) , for instance, find A v i r = 
560 at Zcoii = for the w = —1/3 model instead of 400. The 
general form for the expression of the potential energy in 
the virial theorem yields a weaker dependence of A v i r on w, 
making i t more difficult to constrain w using weak lensing 
(see also lBartelmann. Perrotta fc Baccigalupill2002F) . 



5 ABUNDANCE OF RICH CLUSTERS 

Let us now calculate the comoving number density of viri- 
alised clusters in a certain mass range as a function of red- 
shift: 



n(M,z,w)dM = 



po dv(M, z, w) 



f{v)dM 



M dM 

where po is the background density at redshift zero, and 
8 c (z, w) 



(29) 



v(M,z,w) = 



(30) 



a(M,z,w) 

where S c is the linear overdensity of a perturbation collaps- 
ing at redshift z and a is the r.m.s. of the mass fluctua- 
tion in spheres containing the mass M. S c depends weakly 
on the cosmological model and is close to the Einstein-de 
Sitter value of 1.68. We have ta ken the fit provided by 
IWeinberg fc Kamionkowskil J2003) to S c as a function of z 
and w. 

The standard Press-Schechter mass function is of the 
form: f(y) — y^exp ( — ^ 2 /2). It provides a good general 
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representation of the observed distribution of clusters. How- 
ever, it is known to predict too many low-mass clusters and 
too few high mass clusters, as well as too few clusters at 
high redshifts. We use the modified Press-Schechter func- 
tion proposed by Sheth & Tormen (1999) that gives substan- 
tially better fits to simulated mass functions l|.Tenkins et alJ 
2001). This formula has been given a physical justification 
in terms of an e llipsoid al model for the collapse of pertur- 
bations l|Sheth. Mo fc Tormenll200ll) . The mass function is 

f{v) = y^0.2709(l + 1.1096^ ' 6 ) exp (-0.707// 2 / 2 ) (31) 

We need to calculate the amplitude of the mass fluctu- 
ations a(M), as well as its derivative with respect to M. In 
a Gaussian density field, the variance depends on the scale 
R: 



k 3 P(k)W 2 (kR)=^. 



,1/3 



(32) 



where R(M) = ( ) '" is the Lagrangian radius 
of a halo of mass M at the present time, W(u) = 
3 (sin(tt) — wcos(w)) /w 3 is the Fourier transform of a spher- 
ical top-hat filter with radius R, and P(k) is the power spec- 
trum of density fluctuations extrapolated to z — according 
to linear theory. Assuming that the baryon density param- 
eter Qb,o <§C J2cdm,o, the CDM power spectrum can be ap- 
proximated by P(k) oc kT 2 (k). As for the transfer function, 
we have used 



T(k) 



ln(l+2.34q) 
2.34 9 

[1 + 3.89<? + (16.1c/) 2 + (5.46g) 3 



-1/4 



with q = fc/[n /i 2 Mpc _1 ] jBardeen et alJll986ft . 

We normalise the power spectrum by setting the value 
of the amplitude of the fluctuation within 8/i _1 Mpc spheres 
in a ACDM model. 

erg can be constrained from observations 
of the cosmic microwa v e b ackground (CMB). 
iDoran. Schwindt fc Wetterichl i200ll) provide a useful 
approximation relating erg in a universe with quintessence 
to the erg in a corresponding ACDM universe with the same 
quantity of d ark energy, and cons istent with the COBE 
normalisation (Bu nn fc WhiteHl997f) . In the case of a dark 
energy component with a constant equation of state, as in 
the models we consider here, the expression is simple and 
precise to about 5% : 



^ (l-n x , Q )- {1+w)/r \l^ (33) 

^8 V T 

where to — r(a = 1) is the conformal age of the universe: 

as can als o be inferred from the obse rved abundance of 
X-ray clusters JWang fc Steinhardtlll993) : 



as 



-7(fi ,u>) 



(34) 



where 7 is a function of fio and w. 

In Fig. [7| we show both normalisations of erg: the 
CMB and the X-ray cluster normalisation. We have taken 
= 0.95. as has to decrease with increasing w to meet the 
observational constraints. The CMB normalisation implies 
a steeper decrease of erg with increasing w. 



The cluster evolution is sensitive both to erg and to w 
and might provide a way to constrain those important pa- 
rameters. However, the degeneracy between w and erg pro- 
duces a near cancellation of the two effects. We have shown 
that, as w increases, structures form earlier and are more 
concentrated. However, because of the accompanying de- 
crease in as, more clusters will form at a later epoch. 

Now we want to calculate a more directly observable 
quantity, namely the number density of clusters above a cer- 
tain mass per unit steradian and per redshift interval: 



d 2 N{z) d 2 V c 



d£ldz 



dfldz 



n(M,z)dM 



where the comoving volume is 

d 2 V c d A (zf(l + z) 2 

= c 

dQ,dz 



H{z) 



(35) 



(36) 



and dA is the angular diameter distance. 

Fig. [8] shows that quantity as a function of redshift for 
the different dark energy models and using both the CMB 
and the cluster erg — w relations. As the threshold mass of the 
cluster catalog increases, the peak of the redshift distribu- 
tion is shifted toward lower redshifts and the number density 
of clusters decreases. The amplitude of the curves is sensitive 
to the erg— w relation, but not the position of the peak in red- 
shift. Note that we have assumed a constant limiting mass 
as a function of redshift. It is possible to carry out a similar 
calculation assuming a dependence of Mum with z, appro- 
priate for a given survey. For instance the expected yield 
of a survey of galaxy clusters using the Sunyaev-Zeldovich 
effect is sensitive t o the mass threshold of detection (e.g., 
iHolderet al]l2000ft . Also, such an approach of course as- 
sumes that the cluster masses are known. In practice, the 
total masses have to be inferred assuming some scaling rela- 
tion with an observable (for instance the temperature of the 
X-ray emitting gas). The evolution of the mass-observable 
relation has to be well-calibrated. In principle, this can be 
achieved in dark energy models with constant w if the sur- 
vey is deep enough so that the clusters counts ca n be binn ed 
both in redshift and in the observable quantity (lHull2003l) . 



6 SUMMARY AND CONCLUSION 

Using the top-hat model and the Sheth- Tormen refinement 
of the Press-Schechter mass function, we have examined the 
effect of the equation of state of dark energy on the evolu- 
tion of spherical overdensities and the abundance of clusters. 
We have considered constant-it; models. We have used a gen- 
eral expression for t he virial theorem in the pres ence of dark 
energy instead of the lWang fc Steinhardtl lll99 8) approxima- 
tion widely used in the literature and valid mainly around 
w = — 1. Using the general expression, one can show that 
a dark energy model with an equation of state w — —1/3 
is equivalent to an open CDM model with the same mat- 
ter de nsity f2o and no dark energy. The lWang fc Steinhardtl 
( 1998) approximation in that regime yields an overestimate 
of the overdensity of the virialised cluster A v i r by as much as 
40%. A flat dark energy model with w = —1/3 differs, how- 
ever, from the corresponding open CDM model because of 
the different curvature, implying a different distance-redshift 



8 Cathy Horellou and Joel Berge 



o 




12 3 

z 

Figure 8. Number density of virialised clusters with mass larger 
than a given mass threshold per redshift interval and square de- 
gree in the three cosmological models with dark energy, using the 
X-ray cluster normalisation (thin lines) or the CMB normalisa- 
tion (thick lines). The solid line shows the ACDM, the dashed line 
the w = —0.8 model and the dotted-dashed line the w = —0.6 
model. 

relation and different number densities of virialised clusters 
per redshift interval and per square degree on the sky. 

The number density of massive clusters is sensitive 
both to the amplitude of the mass fluctuations, as, and 
to the as — w relation. As w increases, as has to decrease 
to match the observations of the cosmic microwave back- 
ground and the observed abundance of X-ray clusters. As w 
increases, clusters form earlier and are more concentrated. 
However, this effect is partly counteracted by the decrease 
of as with w. A recent analysis of the evolution of the X- 
ray temperature of clusters favours a low <rg and a higher 
w: a s = 0.6 6 ± 0.16 with w = -(0.42 ± 0.21) at 68% confi- 
dence level (lHenrvl f2004). A previous similar analysis with 
fewer distant clusters in a flat ACDM mode l had yielded a 
higher value of erg (0.77 ± 0.15, lHenrv||200oT) . Breaking the 



as — w degeneracy will require very sensitive observations 
of the evolution of the cluster abundance together with a 
good knowledge of the mass-temperature relation and the 
limiting mass of the survey. 
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